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Abstract 

We consider singular solutions of the biharmonic NLS. In the L^-critical case, the blowup rate 
is bounded by a quartic-root power law, the solution approaches a self-similar profile, and a finite 
amount of L^-norm, which is no less than the critical power, concentrates into the singularity 
("strong collapse"). In the L^-critical and supercritical cases, we use asymptotic analysis and 
numerical simulations to characterize singular solutions with a peak-type self-similar collapsing 
core. In the critical case, the blowup rate is slightly faster than a quartic-root, and the self-similar 
profile is given by the standing-wave ground-state. In the supercritical case, the blowup rate is 
exactly a quartic-root, and the self-similar profile is a zero-Hamiltonian solution of a nonlinear 
eigenvalue problem. These findings are verified numerically (up to focusing levels of 10^) using 
an adaptive grid method. We also calculate the ground states of the standing-wave equations 
and the critical power for collapse in two and three dimensions. 

1 Introduction 

The focusing nonlinear Schrodinger equation (NLS) 

#i(t,x) + AV' + |V'|'"V' = 0, V'(0,x) =Vo(x) Gii'HiK''), (1) 

where x = (xi, . . . ,Xd) € M*^, and A = Yl'j=i Laplacian, has been the subject of intense 

study, due to its role in various areas of physics, such as nonlinear optics and Bose-Einstein Con- 
densates (BEC). It is well-known that the NLS ([1]) possesses solutions that become singular in a 
finite time ^33j. Of special interest is the critical {a = 2/d) NLS 

i^t(t,x) + AV' + |V'|^/''^ = 0, V'(0,x) =V'o(x) Gi/^(M^), (2) 

which models the collapse of intense laser beams that propagate in a bulk Kerr medium. 

In this study, we consider the focusing biharmonic nonlinear Schrodinger equation (BNLS) 

iV'i(t,x)-A2V; + |V'|'"V' = 0, ^{0,^)=M^)(^H\R''), (3) 
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where is the biharmonic operator. Equation ([3]) admits waveguide (standing-wave) solutions of 
the form ^/^(t,x) = A^/°"e*'^''*i?(Ax), where R satisfies the " standing- wave" equation 

- A'^R{^) - R+\R\^''R = 0. (4) 

The BNLS ([3]) is cahed ''L^-criticar , or simply "critical" if ad = 4. In this case, the L norm 
{^^powef^) is conserved under the BNLS dilation symmetry ip{t,x) i-^ L~'^/"il: {t/ L'^ ,'Si./ . The 
critical BNLS can be rewritten as 

iV't(i,x)-AV + |V'l'/''V' = 0, V'(0,x) =V'o(x) eF2(M'^). (5) 

Correspondingly, the BNLS with < 4 is called subcritical, and the BNLS with fifi > 4 is called 
supercritical. This is analogous to the NLS, where the critical case is ad = 2. 

In [3], Ben-Artzi, Koch and Saut proved that the BNLS ([3]) is locally well-posed in H^, when a 
is in the H'^ -subcritical regime 

< cT d < 4, 

0<cj<^d>4. 

Global existence and scattering of BNLS solutions in the if ^-critical case o" = 4/((i— 4) were studied 
by Miao, Xu and Zhao [25] and by Pausader 129]. The latter work also showed well-posedness for 
small data. The iJ^-critical defocusing BNLS was studied by Miao, Xu and Zhao |24j and by 
Pausader [271 [28]. 

The above studies focused on non-singular solutions. In this work, we study singular solutions 
of the BNLS in H^, i.e., solutions that exist in H'^iW^) over some finite time interval t G [0,Tc), 
but for which lim H^'llj/z = oo. The first study of singular BNLS solutions was done by Fibich, 

t >Tc 

Ilan and Papanicolau |12j, who proved the following results: 

Theorem 1. Solutions of the subcritical (ad < 4J focusing BNLS ([3]) exist globally. 

Theorem 2. Let ||'(/'o|l2 < ^cr, where Per = II-RII27 '^^'^ ^ ground state of ^ with a = A/d. 

Then, the solution of the critical focusing BNLS ([5]) exists globally. 

The simulations in [12j suggested that there exist singular solutions for ad = A and ad > 4, and 
that these singularities are of the blowup type, namely, the solution becomes infinitely localized. 
However, in contradistinction with NLS theory, there is currently no rigorous proof that solutions 
of the BNLS can become singular in either the critical or the supercritical case. 

To the best of our knowledge, the only work, apart from [12], which considered singular solutions 
of the BNLS is by Chae, Hong and Lee [9], who proved that if singular solutions of the critical 
BNLS exist, then they have a power-concentration property. See Section [5.31 for more details. 

1.1 Summary of results 

In this work, we consider singular solutions of the focusing BNLS in the i?^-subcritical regime ([6]). 
Our purpose is to characterize these singular solutions: Their profile, blowup rate, power concen- 
tration, et cetera. 
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In some cases, we assume radial symmetry, i.e., that ■ip{t,x) = ipltjr) where r = ||x| 



+ • • • + x^. In these cases, equation ([3]) reduces to 

iV't(t,r)- A2v^ + |V^|2"V = 0, V(0,r) = V'o(r-), (7) 

where 

is the radial biharmonic operator. Specifically, the critical BNLS ([5]) reduces to 

iMt,r)-A^^ij+\i;f^'^i; = 0, V(0,r) = Vo(r). (9) 

The paper is organized as follows: In Section [21 we use Noether Theorem to derive conservation 
laws for the BNLS. We recall that in the critical NLS, the conservation law which follows from 
invariance of the action integral under dilation leads to the NLS "Variance Identity" , which can be 
used to prove the existence of singular solutions. In Section [2. II we use a similar procedure to derive 
the "Variance Identity" for the critical BNLS, and then generalize it to the supercritical BNLS. 
However, since it is not clear that the "BNLS variance" is positive definite, this identity does not 
lead to a proof of the existence of singular solutions. 

The ground states of the BNLS standing-wave equation ^ were previously computed only in 
the one-dimensional case |12j . since they were computed using a shooting method, which cannot 
be easily generalized to multi-dimensions. In Section [3] we use the spectral renormalization method 
to compute the ground-states of the critical BNLS ([5]) for one, two and three dimensions. The 
calculated ground-states provide the first numerical estimate of the critical power for collapse 
Per = 1 1-^1 1 2 i'^ the two-dimensional and three-dimensional cases, see Sectional Direct simulations 
of the critical BNLS suggest that the constant Per in Theorem [2] is optimal. 

In Section [5] we use rigorous analysis to study the critical BNLS ([5|). The blowup rate is shown 
to be lower-bounded by a quartic root, i.e., || AV'||2 < C {Tc — t)^^'^ . The corresponding bound for 
the critical NLS is a square root, i.e., ||VV'||2 ^ < C (Tc — t)^''^ . We then prove that singular solutions 
converge to a self-similar profile strongly in L^^^'^, for any a in the iif^-subcritical regime ([6]). This 
implies that the singular solutions have the power-concentration property, whereby the amount of 
power that enters the singularity point is at least Per. These rigorous results mirror those of the 
critical NLS. 

Let us denote the location of the maximal amplitude of a radially-symmetric solution by 

?'max(i) = argmax|^/;|. 

r 

Singular solutions are called ^^peak-type'^ when rmax(i) = for < t < Tc, and ^^ring-type" 
when rinax(i) > for < t < Tc. In this work, we use asymptotic analysis and numerics to find 
and characterize peak-type singular solutions of the BNLS equation. Ring-type singular solutions 
of the BNLS will be studied elsewhere [2i 13j. 

In Section [6] we use asymptotics and numerics to show that peak-type singular solutions of the 
critical BNLS collapse with the quasi self-similar profile 
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where the self-similar profile is the ground state R{p)- The blowup rate is shown to be slightly 
faster than the quartic-root bound. This is analogous to the critical NLS, where the blowup rate 
of peak-type solutions is slightly faster than the square-root bound, due to the loglog correction 
(the "loglog law"). It is an open question whether the correction to the BNLS blowup rate is also 
a loglog term. 

In Section [7] we consider peak- type singular solutions of the supercritical BNLS. In this case, 
asymptotics and numerics show that singular solutions are of the quasi self-similar form 



the blowup rate of L{t) is exactly a quartic-root, and the self-similar profile B{p) is different from the 
ground-state. Rather, as in the supercritical NLS, the self-similar profile B is the zero-Hamiltonian 
solution of a nonlinear eigenvalue problem. Although B{p) is not in L?, it can be the self-similar 
profile of a collapsing solution, since the collapsing solution is "only" quasi self-similar. 

Section [8] presents the numerical methods. The computations of singular BNLS solutions that 
focus by factors of 10^ and more, necessitated the usage of adaptive grids. We develop a modi- 
fied version of the Static Grid Redistribution method |31i ilOj . which is more convenient for the 
biharmonic problem. Calculating the BNLS standing waves in multi-dimensions is done using the 
Spectral Renormalization Method. 

1.2 Discussion 

In this study, we use rigorous theory, asymptotic theory and numerics to analyze singular solutions 
of the BNLS. All the results presented in this work mirror those of the NLS, "up to a change 
by a factor of 2" in the blowup rate (1/2 — > 1/4), in the critical value of a (2/d — > l/d), et 
cetera. However, several key features of NLS theory are still missing from BNLS theory. First, the 
"Variance Identity" for the BNLS cannot be used to prove that singular solutions exist. Second, 
the critical NLS is invariant under the pseudo-conformal ("lens-transformation") symmetry which 
can also be used to construct explicit singular solutions. At this time, it is unknown whether an 
analogous identity for the critical BNLS exists. Third, in critical NLS theory, the self-similar profile 
is known to possess a quadratic radial phase term, i.e.. 



This term represents the focusing of the solution towards r = 0, and plays a key role in the rigorous 
and asymptotic theory of the critical NLS. At this time, we do not know the analogous radial phase 
term for the critical BNLS. 

Finally, we note that a similar "up to a factor of 2" connection exists between singular solutions 
of the nonlinear heat equation, see [17], and the biharmonic nonlinear heat equation, see [7j. For 
example, the norm of singular solutions blows up as (Tc — t)~^^'^ for the nonlinear heat equation, 
and as (Tc — t)~^^^ for the biharmonic heat equation. The "similarity up to a factor of 2", however, 
is not perfect. For example, the self-similar spatial variable is r/yX^c^^OlosC^c^^ fo'^ the 




lim L(t) = 0, 



(10) 
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nonlinear heat equation, and r/\/Tc — t for the biharmonic nonlinear heat equation. Another 
difference between the equations is that singular solutions are asymptotically self-similar for the 
nonlinear heat equation, and truly self-similar for the biharmonic heat equation. In contrast, the 
NLS possesses self-similar singular solutions, whereas for the BNLS it is unknown whether singular 
solutions are truly, or only asymptotically, self-similar. 



2 Invariance 

The BNLS ([3]) is the Euler-Lagrange equation of the action integral 

5 = y Cdxdt, 

where C is the Lagrangian density 

2 i + cr 

Therefore, the conserved quantities of the BNLS can be found using Noether theorem, see Ap- 
pendix|El As in the standard NLS, invariance of the action integral under phase-multiplications ■0 ^ 
e^^ip implies conservation of the "power" ( norm ), i.e,. 

p(t)^p(o), Pit) = \m)\\l 

Similarly, invariance under temporal translations t ^ t + 5t implies conservation of the Hamiltonian 

H{t) ^ H{0), H[m] = IIA^II^ - ^ \mZXi) ' (12) 

and invariance under spatial translations x x + 5x implies conservation of the linear momentum, 
i.e., 

P(t) =P(0), P(t) = J Im{'0*VV'}dx. 

In the critical case a • d = 4, the action integral is also invariant under the dilation transforma- 
tion 'ip{t,x) X^^'^^lJ{X'^t, Xx). The corresponding conserved quantity is 

J(t) = J(0), J{t) = j x-lm {-0* VV'} dx + UH. (13) 
2.1 Towards a variance identity 

We recall that the action integral of the critical NLS ([2]) is invariant under the dilation transfor- 
mation 

V'NLs(t,x) ^ A^/2^NLs(A^t,Ax). 

The corresponding conserved quantity is 

JNLs(i) = ^NLS(O), JnLS = / x • Im{?/'*VV'} dx - 2tFNLS • (14) 
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In addition, the integral term J x • Im {?/^*VV'} dx is the time-derivative of the variance, i.e., 
-yNLS = 4 J x-Im{V^*V^}dx, V^Ls{t) = J |x|2|V;|'dx. 



di 



Therefore, it follows that 

In the supercritical NLS, the second derivative of the variance is not related to this conservation 
law. Nevertheless, direct differentiation shows that 

^^NLS = 4 J x-Im{^*VV}cZx, 

1 1 x.im#*vv'}dx=2/7NLs-^^iivii2S::;] . 

Therefore 

^FNLS = 8i/NLS-4^M^S:t;]. (15) 



Since Vnls ^ 0, the variance identity (jlSp shows that solutions of the critical and supercritical 
NLS, whose Hamiltonian is negative, become singular in a finite time [341 . 

We next extend the analogy between (I13p and (I14p to the non-critical case. In the case of the 
BNLS, direct differentiation shows that 

4 / x-Im{V^*VV'}rfx = 4g- ll^ll^!"!!!- 
dt J 2(0- + 1) ii^ii2(<7+i) 

Therefore, if we define 



^BNLs(i) = 143NLS(0) + ^ ^(^J X ■ lui {^{3 , x)V ^P} dx^ ds , (16) 



where I^nls(O) is a positive constant, we get the BNLS variance identity 

4w = 4i^ --^(-^^ 



di2''"^^^^ 2(CJ+1) "'^"2('^+l)' 

In order to use this identity to prove singularity formation, however, one must show that I^nlS; 
as defined by (jl6p . has to remain positive. Direct integration by parts gives that 

y"x.Im{V^*V^}dx = ^^^^y |x|4.ImW*AVV}(ix + ^^^i^ Ix^^prfx 

While the second term on the RHS is a temporal derivative of a positive-definite quantity, the first 
term on the RHS is not. Therefore, it remains an open question whether Vbnls has to be positive. 
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3 Numerical calculation of standing waves 

The BNLS equation ([3]) admits the standing-wave solutions 

i^{t,x.) = A2/-e^^'*i?(Ax). 
The equation for the standing-wave profile is 

- A^R{x) - R+\R\^^R = 0, (17) 
where R G H^. For example, in one dimension, equation (|17p is given by 

- Rxxxxix) -R+ m^^R = 0, (18) 

and in two dimensions by 

— {Rxxxx{x, y) + 2Rxxyy + Ryyyy) — R + \R\'^'^R = 0. (19) 

If we impose radial symmetry, eq. pT|) reduces to 

-A^i?(r)-i?+|i?|2'^i? = 0, (20a) 

where is given in ([8]). At r = 0, all the odd derivatives vanish, and so the solution of ([8]) is 
subject to the boundary conditions 

R'{0) = R"'{0) = R{oo) = R'{oo) = 0. (20b) 

The solution of eq. (I20p was computed numerically in the one-dimensional case in jl2] as follows. 
In the ID case, eq. (j20p can be integrated once, yielding an explicit relation between -R(O) and R"{0). 
This parameter-reduction enables the usage of a one-parameter shooting approach. Unfortunately, 
such a parameter reduction is not possible in higher dimensions. Therefore, in multi-dimensions we 
compute the ground-states using the Spectral Renormalization method (SRM), which was intro- 
duced by Petviashvili in [30], and more recently by Albowitz and Musslimani in [Ij. See Section [8.21 
for further details. 

In Figure [1] we display the results in the critical case, i.e., the ground states of 

- A'^R- R + \Rf/'^R = 0. (21) 

Figure [lA] displays the ground-state of eq. in the critical ID calculated by the SRM. 

The solution is in excellent agreement with the solution computed in [12] using the shooting method. 
Figure [TB] and Figure [TC] display the ground state in the critical 2D and 3D cases. We note that, 
while the SRM method that we use does not enforce radial symmetry, the calculated ground states 
for d = 2 and d = 2> are radially symmetric (data not shown) . As noted in [12] , the ground-states of 
the BNLS are non-monotonic in r and change their sign, in contradistinction with the ground-states 
of the NLS which are monotonically-decreasing and strictly positive. 
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5 10 15 5 10 15 5 10 15 

X r 1" 



(A) d = 1 (o- = 4) (B) d = 2 (cr = 2) (C) d = 3 (a = 4/3) 

Figure 1: Numerical solution of eq. (|2ip . using the SRM method (solid line). Red dotted line 
in (A) is the solution computed in [12] using the shooting method. 



4 Critical power for collapse 

Theorem [2] shows that the critical power for collapse in the critical BNLS ^ is Per = ||P"^ 



2' 



when R is the ground-state of (j2ip . The computation of R, see Section [3l allows for the numerical 
calculation of the critical power Per- The case d = 1 was found in [12j to be 



POO 

P^^{d = l)= |P(x)|^dx « 2.9868. 

J x=—oo 



Using the calculated ground state in the two-dimensional case, see Figure \TB\ we now calculate the 
critical power in the two-dimensional case, giving 



p^^{d = 2)= // \R{x,y)\'^dxdy ^13. U3. 

J J x,y=—oo 

Similarly, using the calculated ground state in the 3D case, see Figure IICI gives 



P„(d = 3)= /// \R{x,y, z)\'^dxdydz ^ A4.88. 

J J J x,y,z=—oo 

We now ask whether Theorem [2] is sharp, in the sense that for any e > 0, there exists an initial 
condition t/^Q G such that HV'olli — (l + £)-fcr and the corresponding solution of the critical BNLS 
becomes singular. As noted, at present there is no proof that solutions of the BNLS can become 
singular. Therefore, in particular, it is unknown whether Theorem [2] is indeed sharp. Hence, we 
will explore this issue numerically. 

We recall that in the critical NLS, the necessary condition for collapse HV'lli — -fcr^'^^ := 
||^NLS||2 -g g];jg^].p^ since for any e > the initial condition ipQ = {1 + e)R^^^ becomes singular in 
a finite time [36]. Therefore, we now check numerically whether for < e ^ 1, the solution of the 
critical BNLS with the initial condition = (1 +e)R{x) becomes singular. To do this, we solve 
the one-dimensional and two-dimensional critical BNLS equations with the initial condition = 
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(A) d = l,cr = 4 



(B) d = 2,a = 2 



Figure 2: Collapsing solutions of the critical BNLS ([9]) "with the perturbed ground-state initial 
condition ipo{x) = 1.001 • i?(x). 



1.001 • R{x) , see FigO In both cases, the solution appears to blo"w up, suggesting that Theorem[2] 
is sharp. See Section [8.11 for the numerical methodology used for solving the BNLS. 

We also note that in the critical NLS, if the initial condition is different from the ground state, 
the po"W"er input required for collapse is strictly larger than Pcr'^^^i see [TT]. This the case also 
in the BNLS. For example, in the one-dimensional critical BNLS "with a Gaussian shaped initial 

2 

condition ipQ = C ■ , the input po"wer required for collapse is strictly larger than 1.003 • Pa, 
see |12] . In the t"wo-dimensional critical BNLS "with a Gaussian shaped initial condition, the input 
power required is strictly larger than 1.001 • Per, see [20j . 

5 Blowup rate, blowup profile, and power concentration (critical 
case) 

5.1 Lower-bound for the blowup rate 

In [8j, Cazanave and Weissler proved that the blowup rate for singular solutions of the critical 
NLS (l2|) is not slower than a square-root, i.e., that ||V'0||2 ^ K{Tc — t)~^/'^ . The analogous result 
for the critical BNLS is as follows: 

Theorem 3. Let ijj he a solution of the critical BNLS (0j that becomes singular at t = T^ < oo, 

— 1/2 

and let l{t) = ||A'0||2 . Then, 3K = i^(||'0o|l2) > such that 

l(t) <KiT^-tf/'^, 0<t<Tc. 

Proof. We follow the proof given by Merle [21] for the critical NLS. For a fixed t, < t < Tc, lei 
us define 

V'i(s,x) = f/^ilj{t + S-f, X- /). 
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Then, is defined for t + I s < T^ <^=^ s < Sc = I (t) ■ [Tc — t), and satisfies the BNLS equation 



iS.^i + A^V^i + lV'il^/'^Vi = 0. 



Since 



||A^i||2 = /4 ||AV'(t + s- Z*, X- /)| 

this imphes that hms^^^ ||A'(/'i||2 = oo, i.e., that V'i(s) becomes singular as s — > 5c- In addition, 

\\AMs = 0,x)||2 = ||AV'(t,x)||2 = 1. (22a) 

Prom the definition of ipi and power conservation it follows that 

||Vl(. = 0,x)||^ = ||V(^,x)||^ = ||^o(x)||^ (22b) 

Using equations ()22ap and ()22bp and the Cauchy-Schwartz inequality gives 

\\VMs = 0,x)||2 < = 0,x)||2 • \\AMs = 0,x)||2 = ||Vo(x)||2 . (22c) 

Together, the three formulae (|22|) imply that for any fixed t £ [0,Tc), 

||Vi(s = 0,x)||^2 < llV^oll^ + llVolls + l- (23) 

In other words, for each t, the initial norm of ipi is bounded by a function of ||V'o|l2- Specifically, 
this bound is independent of t. Prom the local existence theory ipi exists in s G [0,SM{t)], 
where Sm = Sm = ^,^)\\h^) ■ Therefore, it follows from (p3]) that Sm depends on HV'oll) but 

is independent of t. 

Since ipi blows up at 5c we have that 

Sm <Sc{t) = l-Ht)-{T,-t), 

from which the result follows- □ 

5.2 Convergence to a quasi self-similar blowup profile 

In [37], Weinstein showed that the collapsing core of all singular solutions of the critical NLS 
approaches a self-similar profile- We now prove the analogous result for the critical BNLS: 

Theorem 4. Let d > 2 and let ip{t, r) he a solution of the radially- symmetric critical BNLS ^ 

with initial conditions ipo^r) £ ^^radiai' ^^'^^ becomes singular at t = Tc < oo. Let l{t) = \\Aip\\2 
and let 

Smt,r) = f'\t)i,{t, l{t)r). 

Then, for any sequence t!^ Tc there is a subsequence tk such that S{ilj){tk,r) ^'(r) strongly 
in L'^, for all q such tha 



|2<g<oo 2<d<4, 

1 1 2 1 1 2 

In addition, \\^\\2 > where R is the ground state of equation (^j). 



2<g<^ 4<d. ^^^^ 



^ In fact, q — 2{a + 1), where a is in the //^-subcritical regime ((6)1. 
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Proof. Let ^ T^ and define 

0fc(r) = Smtk^r) = f'\tMtk,l{tu)r). 
From the definition of (pk it follows that 

UkWl = UoWl , W^Ml = = 1, H[^k] = l^H[^Pit,)]. 

Therefore, using Cauchy-Schwartz, 

||V0fc||^ < ||(^fc||2- ||A,/.fc||2 = ll^olla- 

Since ||</'A:||j:/2 is bounded, it follows that there exists a subsequence of (pk which converges weakly 
in to a function ^ G -f^^adiai- From the Compactness Lemma [131 see Appendix [A) it follows 
that (pk ^ ^ strongly in L"?, for all q given by ()24p . 



Next, we prove that H[ii] < 0. Since (pk^^, it follows that A(j)k^A'$, and so ||AvI'||2 < 
limfc_»oo ||A</)A;||2 = 1- Additionally, since (pk ^ ^ for some g = 2 + 2s > 2, we have that ||^'||2 = 
limfc_^oo \\Pk\\2, and so 

H[<if] < lim H[pk] = lim I'^Hltpo] = 0. 

In addition, since 

= lim H[Pk] = lim (l - mXV^^ , 

it follows that limfc^oo ll</'fc|l2(CT+i) > 0, so ^' 7^ 0. Therefore, Corollary [TSl see Appendix [Q implies 
that ||^'||2 > \\R\\l. □ 

5.3 Power Concentration 

Solutions of critical NLS have the power concentration pro|)erty, whereby the amount of power that 
collapses into the singularity is at least Pcr^'^^ = ||-^^'^^||2' [371 123]. Iii what follows, we prove 
the analogous results for the critical BNLS. 

Corollary 5. Let d>2, and let ip{t,r) he a solution of the radially- symmetric critical BNLS ^ 
that becomes singular at t = Tc < 00. Then, Ve > 0, 

liminf ||VX*,r)llL2(r<.) > ^cr, 

where Per = II-RII2 • 

Proof. The result shall follow directly from Corollary [6l □ 

The following Corollary shows that the rate of power-concentration is not slower than the 
blowup rate l{t). The NLS analogue is due to Tsutsumi [35j and Weinstein |37j . 
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Corollary 6. Let d > 2, let ip{t, r) he a solution of the radially- symmetric critical BNLS (0j that 

— 1/2 

becomes singular at t = Tc < oo, and let l{t) = ||A^||2 . Then, 

1. For any monotonically- decreasing function a{t) : [0, Tc) such that 

lim a{t) = 0, and lim //a = 0, 

we have that 

liminf ||V'(t,r)||2 2(^^^(^^^ > P„. 

2. For any e > 0, 3K > such that 

limm{\\'tP{t,r)\\l2(^^^Kl{t)) > (1 " e)^cr. 

Proof. See Appendix [Bl □ 

Since the second part of Corollary [6] is true for all e > and since Kl{t) 0, Corollary [5] follows. 

The next Corollary shows that the power-concentration rate has a quartic-root upper bound. 
The analogue in NLS theory, which is a square-root upper bound, was proved in 



Corollary 7. Let d > 2 and let ilj{t,r) he a solution of the radially- symmetric critical BNLS ^ 
that hecomes singular at t = Tc < oo. Then, 

1. For any monotonically- decreasing a : [0,Tc) M"*" such that 

- ty/^ 

lim a(t) = 0, and lim — = 0, 

t-^T, ^ ' t^n a{t) 

we have that 

limmf ||?/;(f,r)||2 2(^^^(^j^ > Per- 

2. For any e > 0, > such that 

liminf ||V'(t,r)||^2, ^^(y^_^)i/4) > (1 - e)Pcr. 

Proof. Theorem [3] implies that 1(f) < K{Tc — i)^/^. Therefore, the result follows immediately from 
Corollary [6l □ 

Recently, Chae, Hong and Lee |9j used the harmonic analysis method of Bourgain ^ to prove 
that singular solutions of the critical BNLS for d>2 have the power-concentration property 

lim supj|^(t,x)||2^^|^_^^|^^^^_^^,/,^ >C7, 



Xo 



where C is a positive constant. This result is more general than Corollaries I5|6|7l in that it does not 
assume radial symmetry. The proof given here, however, is considerably simpler. More importantly, 
it shows that C = Per- 
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6 Peak-type singular solutions of the critical BNLS 

In this, we consider radially-symmetric singular solutions that are ^^peak-type" , i.e., for which r^ax{t) = 
for < t < Tc, where rmax(i) = argmax|^| is the location of the maximal amplitude. 

r 

6.1 The critical NLS - review 

The critical NLS admits singular solutions that collapse with the universal ippjNLs profile, i.e., V ~ 
V'ijNLs, where 

1 -?NLS^^\ ir+j|ir2 /"* ds 



The self-similar profile R^^^ is the ground-state solution of 

-R + AR + lRf" R = 0. 
The blowup rate of L{t) is given by the loglog [H I [T8 | [T9| [22] 

f 2tt{T, -t) \^ 

Lit) ~ r -= -r, , t^Tc. (26) 

Viog|iog(rc-t)|y 

Since the blowup rate ()26p is slightly faster than a square root lim LLi = lim — (^^)^ = 0. 

Therefore, the phase term jj-r"^ = -^/O^ in ([25]) vanishes as t ^ Tc- Hence, the blowup profile 
reduces to 

6.2 Informal analysis 

We now look for the "corresponding" peak- type singular solutions of the critical BNLS (0). The- 
orem m suggests that the collapsing core of the singular solution approaches a self-similar form, 
i.e., 

2Pit,r) r^tPB{t,r), 0<r<p,-L{t), 

where 

^B{t,r) = j^^B{p)e'^^'\ p=^, (28) 
and pc = 0(1). Substituting ([28]) into ^ and requiring that [ipt] ~ [AV'] ~ [IV'I'^'^^V'] suggests that 

Let us consider the self-similar profile B{p). In the singular region r = 0{L) we have that 
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and 



i^t - - [^B - L,L' (^S + pB,^ } . 

Hence, B{p) satisfies 

- Bip) - A^B + \Bf' = i (^lim LtL^^ (^^B + pB^^ . (29) 

Theorem [3] shows that that blowup rate of L{t) is lower-bounded by a quartic-root. In the 
critical NLS the blowup rate of peak-type solutions is slightly faster than the analogous square- 
root rate, due to the loglog correction. Hence, we expect that the blowup rate of peak-type critical 
BNLS solutions is slightly faster than a quartic root, i.e., 

^<*> 0. (30) 



In that case, lim LfL^ = lim — (L'^)^ = 0, and (1291) reduces to the standing- wave equation (1211) . 

t^Tc t^Tc 4 ^ ^ * ^ 

Since the ground-states of ([2T]) attain their maximal amplitudes at p = 0, see Section [3l these are 
peak-type solutions. 

The above informal analysis thus leads to the following Conjecture: 

Conjecture 8. The critical BNLS admits peak-type singular solutions such that: 

1. The collapsing core approaches the self-similar profile 

ijit, r) ~ ijR{t, r), < r < pe • Ht), (31a) 

where 

and R is the ground-state of equation (I21|) . 

2. The blowup rate of L{t) is slightly faster than a quartic-root, i.e., 

hm-^=|° ^="\ (32) 
t^T,{T^-t)P |cx) p>l/4 ^ ' 

In Section [6T3l we provide numerical evidence in support of Conjecture [H) 

6.3 Simulations 

The one-dimensional critical BNLS 

i^t(t, j;) - Vxxxx + IV'I^V' = (33) 
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le-07 2e-07 
t 



(A) d = 1 (B) d = 2 

Figure 3: Maximal amplitude of peak-type singular solutions of the critical BNLS Q. 



was solved with the Gaussian initial condition ipo{x) = A\e.~^ with A\ ~ 1.618, whose power 
is 11^0 II 2 — 1-1 ■ -fcr('^ = 1)- The maximal amplitude of the solution HV'lloo ^ ^ function of time 
is plotted in Fig I3A[ The amplitude increases abruptly by a factor of 10^ around T^ ~ 0.0499, 
suggesting that the solution becomes singular in a finite time. 

The simulation was repeated for the radially-symmetric two-dimensional critical BNLS 



1 



1 



-A 



Iprrrr + 4' = 0, 



(34) 



with the Gaussian initial condition ipo{r) = A2e~^ with A2 ~ 3.034, whose power is HV'olli 
l.lPcr{d = 2) The amplitude increases abruptly by a factor of 10^. around Tc ~ 0.0606, see Fig[ 
again suggesting that the solution becomes singular in a finite time. 

We next consider the self-similar profile of the collapsing solutions from Figure [3j In order to 
verify that it is given by (i3T]l . we rescale the solutions as 



V'rescaled(t,p) = L^/^ (t) ^^(t , r = /> • L) , L{t) 



-0-/2 



(35) 



with 2/a = d/2. The rescaled solutions at focusing levels of L = 10^^ and L = 10~^ are indistin- 
guishable, see FigHl indicating that the collapsing core is indeed self-similar according to (j3ip . A 
predicted, the self-similar profile is very close to the ground-state R in the core region < p < 4. 
We next compute the blowup rate p, defined by the relation 

L ~ k{t, - typ. 

To do that, we perform a least-squares fit of log(L) with log(rc — t), see Figure El obtaining a 
value of p ~ 0.2516 for both d = 1 and d = 2. This value of p is slightly above 1/4, implying that 
the quartic-root lower-bound given by Theorem [3] is close to the actual blowup-rate of peak-type 
singular solutions. 

Next, we provide two indications that the blowup rate is faster than 1/4. First, if the blowup 
rate is exactly 1/4, then lim^^T^^ L^Lt should be finite and strictly negative. However, up to focusing 
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x/L r/L 

(A) d = 1 (B) d = 2 

Figure 4: The solutions of Figure [3l rescaled according to (|35p . at focusing levels L{t) = 10~^ 
(blue dotted line) and L{t) = 10~^ (black solid line). Red dashed line is the rescaled ground- 
state The three curves are indistinguishable for < r/L < 4. 




(A) d = 1 (B) d = 



Figure 5: L{t) as a function of (Tc — t) , on a logarithmic scale, for the solutions of Figure [3] 
(circles). A) ID case. Solid line is the fitted curve L = 0.742 • (Tc - t)°-25i6 . b) 2D case. 
Sohd line is the fitted curve L = 0.641 • (Tc - t)0-25i6 . 
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Figure 6: A: L^Lt as a function of 1/L, for the solution of Figure [3Al (black solid line) and 
of Figure [3B] (red dashed line). B: same as (A) for the supercritical cases d = l,a = 6 (black 
solid line) and d = 2,a = 3 (red dashed line). 



level of L = 10~^, L^Lt does not appear to converge to a negative constant, but rather to increase 
slowly towards 0~, see Figure [6Al Second, according to the informal analysis in Section [6.2^ the 
blowup rate is faster than a quartic root if and only if the self-similar profile B{p) satisfies the 
standing- wave equation (jl]), which is indeed what we observed numerically in Figure HI 

Remark: In the critical NLS the blowup rate of peak-type solutions is slightly faster than the 
analogous square-root lower-bound, due to the well-known loglog-correction (1260 . Figure [5] shows 
that the blowup rate is slightly faster than a quartic root, and Figure I6AI shows that L'^Lt — > 
very slowly. Together, this suggests that the blowup rate in the critical BNLS is only slightly 
faster than the analogous quartic root. At present, we do not know if the blowup rate of peak-type 
solutions of the critical BNLS is a quartic root with a loglog correction. We note, however, that 
the loglog correction in the critical NLS cannot be determined numerically [13], and can only be 
derived analytically. Therefore, we expect that the determination of the analogous correction to 
the 1 /4 blowup rate of the critical BNLS will also have to be done analytically, and not numerically. 



7 Peak-type singular solutions of the supercritical BNLS 
7.1 The supercritical NLS - review 

In contrast to the extensive theory on singularity formation in the critical NLS, much less is known 
about the supercritical NLS. Numerical simulations and formal calculations (see, e.g., J33, Chapter 
7] and the references therein) suggest that peak-type singular solutions of the supercritical NLS 
collapse with a universal V'Q profile, i.e., 

^ \^o{t,r) 0<r<rJ, 

I V'non-singular (^) '") ^ ^ 
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where 

Note that the singular region r G [0, rd is constant in the coordinate r. Therefore, in the rescaled 
variable p = r/L{t), the singular region p G [0,rc/L{t)] becomes infinite as L{t) — > 0. This is in 
contradistinction with the critical case, wherein the singular region p E [0, pd is constant in the 
rescaled variable p, but shrinks to a point in the original coordinate r. 
The self-similar profile Q is the solution of 

Q"{p) + ^Q' -Q + i'^(-Q + pQ') + \Q\'^Q = o, 

p 2 \a ) (36c) 

Q'(0) = 0, Q(oo) = 0. 

Solutions of ()36cp are complex- valued, and depend on the parameter k and on the initial condi- 
tion (5(0) = Qo- Solutions of (j36cp whose amplitude IQI is monotonically-decreasing in and 
which have a zero Hamiltonian, are called admissible solutions [33]. For each choice of {cr,d), equa- 
tion (j36cp has a unique admissible solution (up to a multiplication by a constant phase e*"). This 
solution is attained for specific real values of k and (5(0), which we denote as 

K = KQ{a,d), Q{0) =Qo{a,d). (37) 

The blowup rate of L{t) is a square-root, i.e.. 



L{t) ~ K^/Tc - t, t ^ Tc. (38) 
Numerical simulations and formal calculations suggest that: 

1. The self-similar profile of singular peak- type solutions of the NLS ([1]) is an admissible solution 
of ()36cp . Since Q{p) attains its maximal amplitude at p = 0, the solution is peak-type. 

2. The constant k of the blowup rate ([38|) is equal to hq[(t, d). Hence, in particular, k is universal 
(i.e., is independent of the initial condition tpo). 

The admissible solution Q{p) satisfies 

|Q(p)| ~C-p-V-, p^oo. 

Thus, Q ^ L2(R). Nevertheless, Q{p) can be the self-similar profile of solutions, since il}{r^t) ~ 
V'Q only for r G [0,rc], see [5j. 

7.2 Informal analysis 

As in the supercritical NLS, we expect that singular peak-type solutions of the supercritical BNLS 
collapse as 

V'(i,n ~ < , . (39) 

I V^non-singular v^) ^ j ^ ^ 
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where ■0B is a self-similar profile, to be determined. As in the supercritical NLS, the singular 
region r G [0, rd, is constant in the coordinate r. Therefore, in the rescaled variable p = r/L(t), the 
singular region p G [0,rc/L{t)] becomes infinite as L(i) — > 0. This is again in contradistinction with 
the critical-BNLS case, where the singular region p G [0, pd is constant in the rescaled variable p, 
but shrinks to a point in the original coordinate r. 

The BNLS ([3]) is invariant under the dilation symmetry ri-^ -£,^1— > -^,■01— > jj/^ip, where L is 
a constant. In the supercritical case ad > 4, this suggests that 

^Bit,r) = j^^B{p)e-^^\ p=L. 

Similar arguments as in Section [6.21 show that r(t) = jl^Q l^{s) ^^' Therefore, as in the supercritical 
NLS, we expect the collapsing part of the solution to approach the self-similar profil^ 

Theorem [3] showed that in the critical case, if L{t) ~ k(Tc — t)^, then p > 1/4. The following 
Lemma extends this result to peak- type solutions of the supercritical BNLS. 

Lemma 9. Let ad > 4, and let be a peak-type singular solution of the BNLS that collapses with 
the ij)B profile (I40p . // L[t) ~ k,{Tc — tY , then p > \. Furthermore, p = 1/4 if and only if the 
self-similar profile B{p) satisfies the equation 

- B{p) +i'^ (^^B + pB'^ - ApB + \B\'^''B = 0, k>0. (41) 

Proof. If ip ~ ipB, then 



and 

L^+s/'^ 1 \a 



~ (^B)t ~ TTT^ i iB - LtL^ ( Ib + pB^ 



Hence, the equation for B is 



-B-i hm LtL^^ {^B + pB'^ - Af^B + {Bl^" B = 0, (42) 

implying that LtL^ should be bounded as t ^ Tc. Since L^L, ~ —pK'^{Tc — t)^^ ^, it follows 
that p > i. If p = 1/4, then L^Lt and equation (gS]) reduces to ([31]). □ 



Note that in the critical case 2/a = d/2, hence the self-similar profile (|4Up reduces to (|3f p. 
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Let us consider the fourth-order nonhnear ODE ()4ip for the self-similar profile B{p). Its solution 
requires four boundary conditions, and the determination of the parameter k. Radial symmetry 
implies that -B'(O) = B"'{0) = 0. Since the solution is invariant up to rescaling B{p) — > X'^^'^ B{Xp), 
one can set, with no loss of generality, B{0) = 1. Therefore, we require two additional constraints 
in order to determine k and Bq. We recall that in the supercritical NLS the "admissible value" 
of K is determined from the requirements that Q has a zero Hamiltonian and a monotonically- 
decreasing amplitude, see Section 17.11 For the BNLS, the following informal argument suggests 
that the zero-Hamiltonian condition should also holds. Indeed, from Hamiltonian conservation 
it follows that H [iI^b] is bounded, because otherwise the non-singular region would also have an 
infinite Hamiltonian. Calculating the Hamiltonian of ^b, we have 



rc/L(t) 

2 ^ I J5|2+2cr \ 

p=0 



\B{pr--^\Br'-]p'-'dp 



From //^-subcriticality, see ([6]), it follows that L ^"^'^ ^ cxd as L ^ 0. Therefore, if H[iPb] 
remains bounded as t ^ then 



f oo 



H[B]= I ( |Ap5|2--^|S|2+2")/-idp = 0. (43) 
Jp=0 \ 1 + cr 



WKB analysis of (HTj) shows that, see Appendix [Pj 

B{p) ~ ciBi{p) + C2B2{p) + C3B3{p) + CiBi{p), p^ oo, 



where 



Bsip) - ^(..-, ^ exp(+.-.//3-.-log(,)) 



P 

3V3 



MP) ~ ^ exp(+i^.//3_,_Lio,(,)) 



and b = . Equation (j4ip therefore has two algebraically-decaying solutions, Bi and B2, an 

exponentially- increasing solution B3, and an exponentially-decreasing solution B4^. Since ad > 4, 
the exponent ^{ad — 1) of B2 is larger that the exponent ^ of hence Bi ^ B2 as p ^ 00. 

Lemma 10. Let B(p) be a zero-Hamiltonian solution of (j4ip . Then, C2 = C3 = and 

B{p) ciBi{p), p^oo. 



Furthermore, B" € L^. 
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Proof. The exponentially increasing solution B3 must vanish identically if the integrals are to 
converge, hence C3 = 0. Next, the p^^^ phase term in B2 implies that 

Hence, the integral 

l|-°2 h P P "-P 

diverges in the //^-subcritical regime a{d — A) < 4, i.e., B2 ^ L^. Since, in addition, 

B'l G L^, Bi£ L2+2a^ ^ L2+2<7^ 

the Hamiltonian can be finite only if C2 = 0, in which case B" G L^. □ 

Lemma [TOl shows that the condition H[B] = imposes the two constraints C2 = and C3 = 0. 
Hence, zero-Hamiltonian solutions of equation (|4ip satisfy the five boundary conditions 

B{0) = 1, B'{0) = 0, B"'{0) = 0, C2 {B'^;, k) = C3 (S^,', k) = 0. 

Therefore they form a discrete set of solutions and of values of k. We conjecture that the additional 
condition of monotonicity of will lead to a unique admissible solution B and a unique value 
of K. 

Corollary 11. Let B{p) he a zero-Hamiltonian solution of ()1T]) . Then, II-BII2 = 00. Neverthe- 
less, lim UB\\L2^r<r,) < 

Proof. Since B{p) ~ ciBi{p), 



00 



\\Bi\\l ~ C / p-^/'^+'^-'dp ~ 
Following the arguments of [5], the profile ipB satisfies 

rrc/L{t) 



00 

= 00. 

p=0 



Jp=0 



dp 



□ 



In summary, we conjecture the following: 
Conjecture 12. 

Let be peak-type singular solution of the supercritical BNLS. Then, 
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Table 1: Position of circles in Figure O 



1. The collapsing core approaches the self- similar profile ipB, i-^-, 

ipit, r) ~ ipBit, r), <r < r^ 

where 

r{t) 



2. The self-similar profile B{p) is the solution of 

-B{p) + 05 + pB'^ - ApB + = 0, 

B{0) = 1, B'{0) = B"'{0) = 0, H[B] = 0, 
where k > and H[B] is the Hamiltonian of B, see (j43p . 

3. In particular, B{p) ^ R{p)- 



-ds. 



(44a) 
(44b) 



(44c) 



4- Equation (I44cp has a unique "admissible solution" with a unique "admissible value" of k 
K{a,d), such that \B{p)\ is monotonically decreasing. Additionally, B{p) 



CO. 



5. The admissible solution is the self- similar profile B oftps, see (I44b[) . 

6. The blowup rate of singular peak-type solutions is exactly a quartic root, i.e., 



L{t) ~ K VTc - t, K > 0. 



(44d) 



7. The coefficient k of the blowup rate of L{t) is equal to the value of k of the admissible solu- 
tion B, i.e., 

— K[a, d). 



K := 



lim 



In particular, k is universal (i.e., it does not depend on the initial condition). 
In Section [7.31 we provide numerical evidence in support of Conjecture I12[ 
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(A)d=l,cr = 6 (B)d = 2,a = 3 

Figure 7: Maximal amplitude of peak- type singular solutions of the supercritical BNLS. 



T 2/0, I 
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(A) d = l,a = 6 
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T 2/a, I 




(B) d = 2,0- = 3 



Figure 8: The solutions of FigureO rescaled according to ([35]) . at the focusing levels 1/L = 10^ 
(blue solid line) and 1/L = 10^ (black dashed line). The magenta dotted line is the rescaled 
ground-state R. 
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(A)d=l,CT = 6 (B)d = 2,(j = 3 

Figure 9: The solutions of Figure[71 rescaled according to (j35p . at focusing the level 1/L = 10^ 
(circles). Solid lines are the fitted curves y = 0.63 • {r/L)~^'^^ (left) and y = 0.85 • (r/L)""'^^ 
(right). 



7.3 Simulations 

The radially-symmetric BNLS ([7]) was solved in the supercritical case d = l,a = 6 with the 
initial condition ipo{x) = 1.6e~^ , and in the supercritical case d = 2, a = 3 with the initial 
condition ipo^r) = 3e~^ . In both cases, the solutions blowup at a finite time, see Figure [71 

To check whether the solutions collapse with the self-similar profile (j44p . the solution was 
rescaled according to (135]) . The rescaled solutions at the focusing levels L = 10~^ and L = 10~^ 
are indistinguishable in both the one-dimensional case (Figure l8Ap and the two-dimensional case 
(Figure ISBp , providing numerical support that the solution collapses with the ip b profile (j44p . As 
predicted, the self-similar profile is different than the ground-state R. Indeed, Figure [9] shows that 
as /3 — > oo, the self-similar profile of ip decays as /0~^/°", which is in agreement with the decay rate 
of B,{p). 

We next verify that the solution converges to the asymptotic profile for r G [0, Tc], i.e., for p S 
[0,rc/ L{t)]. To do this, we plot in Figure fTOl the rescaled solution at focusing levels of 1/L = 
10,100,1000,10000, as a function of log(r/L). The curves are indistinguishable at r/L = 0(1), 
but bifurcate at increasing values of r/L. These "bifurcations positions" are marked by circles in 
Figure [TOl and their r/L values are listed in Table [H The "bifurcation positions" are linear in 1/L, 
indicating that the region where ip ~ V'B is indeed p £ [0,rc/L{t)], which corresponds to r G [0, Tc]. 

To compute the blowup rate p, we performed a least-squares fit of log(L) with log(Tc — t), 
see Figure [TT] The resulting values are p « 0.2502 in the d = l,a = 6 case and p « 0.2504 in 
the d = 2, o" = 3 case. Next, we provide two indications that the blowup rate is exactly 1/4, i.e., 
that 

L(t) ~ Ky^Tc - t, K>0. 

First, if the blowup rate is exactly a quartic root, then L^Lt < 0. Indeed, Figure [6B] shows 



7 PEAK-TYPE SINGULAR SOLUTIONS OF THE SUPERCRITICAL BNLS 



25 




(A) d = l,cr = 6 




(B) d = 2,a = 3 

Figure 10: Convergence to a self-similar profile. The solutions of Figure [71 rescaled according 
to (f35]) . as a function of log(r/L), at the focusing levels L = 10^^ (dashed blue line), L = 10^^ 
(dash-doted red line), L = 10^^ (dotted green line), L = (solid black line) and L = W^^ 
(solid magenta line). The circles mark the approximate position where each curve bifurcates 
from the limiting profile, see also Table [TJ 
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(A)d=l,a = 6 (B)d = 2,(j = 3 

Figure 11: L(t) as a function of (Tc — t) , on a logarithmic scale, for the solutions 

of Figure [7] (circles). Solid lines are the fitted curves L = 1.048 • (Tc - (A) 
and L = 0.931 • (T^ - t)0-2504 (b). 



that in the case d = 1, a = 6, L^Lt —0.289, implying that 

K{d = 1, a = 6) f« v^4 • 0.289 « 1.037 . (45) 
In the case d = 2, a = 3, L^Lt -0.171, implying that 

K{d = 2, a = 3) f« v^4 • 0.171 f« 0.909 . (46) 

Since L^Lt converges to a finite, negative constant, this shows that the blowup rate is exactly 1/4. 

Second, according to Lemma [U if L'^Lt the self-similar profile B{p) should not satisfy the 
standing- wave equation . In Figure [8] we saw that the rescaled self-similar BNLS solutions and 
the corresponding ground-states differ considerably (compare with Figure H]) , indicating that the 
blowup rate is exactly 1/4. Therefore, the numerical results again support Conjecture 1121 

Finally, we verified that the value of k in the blowup rate (j44d|l is universal. We solve the BNLS 
in the case d = 1, o" = 6 with the initial condition ^oix) = 2e~^ . In this case, the calculated value 
of K{d = 1, cr = 6) is K = lim y— 4LtL^ ~ 1.037, which is equal, to first 3 significant digits, to the 

t >Tc 

previously obtained value, see (US]), for the initial condition tpo{x) = 1.6e~^ . Similarly, in the case 
d = 2,a = 3, we solve the equation with the initial condition ipo^x) = 3e~^ . The calculated value 
of K{d = 2, cr = 3) is K = lim y— 4LiX^ ~ 0.913, which is equal, to first 2 significant digits, to the 

t 

previously obtained value, see (j46p . for the initial condition "00(3^) = 3e~^ . 
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(A) old method, L = lO"*^ 

Figure 12: The grid-spacing Ar^ obtained using the SGR method of [TO] for a peak-type 
singular solution of the BNLS. A) The grid generated the original method of [TO] at focusing 
level of L = 10~^. The Singular and non-singular regions are well-resolved, but the transition 
region Ar^ displays a discontinuity. At this point, the finite difference operator becomes ill- 
conditioned. B) same as (A), after adding the new penalty function w^, at focusing level L = 
10~^^. Even at this much larger focusing level, the transition region is now well resolved. 
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(B) new method, L = 10~^^ 



8 Numerical methods 

8.1 Adaptive mesh construction using the SGR method 

In this study, we computed singular solutions of the BNLS equation dT]). These solutions become 
highly- localized, so that the spatial scale-difference between the singular region r — r^ax = 
and the exterior regions can be as large as 0{1/L) ~ 10^'^. In order to resolve the solution at both 
the singular and non-singular regions, we use an adaptive grid. 

We generate the adaptive grids using the Static Grid Redistribution (SGR) method, which 
was first introduced by Ren and Wang [31j, and later simplified and improved by Gavish and 
Ditkowsky [10]. Using this approach, the solution is allowed to propagate (self-focus) until it 
becomes under-resolved. At this stage, a new grid, with the same number of grid-points, is generated 
using De'Boors 'equidistribution principle', wherein the grid points {rm} are spaced such that a 
certain weight function is equidistributed, i.e., that 

rrm+i 

/ wi [i^{r)] dr = const, 

J r=rm 

see ISH Ho] for details. 

The algorithms of [311 [TO] keeps a recursive set of grids, the coordinates of each is given in the 
reference frame of the previous grid. This enables the simulation to reach high focusing levels, where 
the grid-points position cannot be stored at the physical reference frame due to loss of significant 
digits. In this study, however, we implement a simplified version of the method of [TO], in which 



8 NUMERICAL METHODS 



28 



we dispense with the above hierarchy of grids and store the grid-points position in terms of the 
original reference frame. On this highly non-uniform grid, we use a standard non-uniform finite- 
difference approximation of the radial biharmonic operator ([8]). The approximation is third-order 
accurate, with a seven-point stencil. Thus, we are limited by the standard machine accuracjH to 
focusing levels of no more than L ~ 10^^^. However, the resulting gain in software (and numerical) 
simplicity justifies this limitation. 

The method in allows control of the fraction of grid points that migrate into the singu- 
lar region, preventing under-resolution at the exterior regions. This is done by using a weight- 
function W2, which penalizes large inter-grid distances. However, we found that this numerical 
mechanism, while necessary, is insufficient for our purposes. In order to understand the reason, 
let us consider the grid-point spacings Ar^ = r^+i — Tm- Using the method of [lOj with both wi 
and W2 causes a very sharp bi-partition of the grid points - to those inside the singular region, 
whose spacing is determined by wi and is Arm = (^{L), and to those outside the singular region, 
whose spacing is determined by W2 and is Ar^ = 0(1), see Figure I12AI Inside each of these 
regions, the finite difference approximation we use is well conditioned. However, at the transition 
between these two regions, the finite-difference stencil, seven-points in width, spans grid-spacings 
with 0{1/L) scale-difference — leading to under-resolution which completely violates the validity 
of the finite-difference approximation. 

In order to overcome this limitation, we improve the algorithm of [lOj by adding a third weight 
function 



which penalizes the second-difference A^r^ = Ar^+i — Ar^ operator of the grid locations, allowing 
for a smooth transition between the singular region and the non-singular region, see Fig I12BI 

On the sequence of grids, the equations are solved using a Predictor-Corrector Crank-Nicholson 
scheme, which is second-order in time. 

8.2 The Spectral Renormalization Method 

Here, we describe the adaptation of the Spectral Renormalization method (SRM) to the standing- 
wave BNLS equation (jl7p . Denoting the Fourier transform of R{x) by (k), equation (|17p 
transforms to 



Typically this iterative scheme diverges either to oo or to 0. In order to avoid this problem, we 
renormalize the solution as follows. Multiplying equation (j47p by J^[R]* and integrating over k 

^ At the time of writing, the standard machine accuracy is IEEE 64-bit floating points, with relative machine 
error of 10~^^. Using 128-bit floating points wiU, for aU practical purposes, eliminate the above restriction. 





(47) 



where fc^ = |k|^, leading to the fixed-point iterative scheme 




m = 0, 1, . . . 
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gives the integral relation: 

SL = SR, where SL[R] = j \T[R]\'^ dk, SR[R] = j -^^^J^ [{Rl"^" R] T[R]*dk. (48a) 

We now define R^,i = CmRm such that the integral relation (I48ap is satisfied by R^,i, i.e, that 

SL [Rm+1/2] = CmSL[Rm] = C^^^SR[Rm] = SR [Rm+1/2] ) 
1 

leading to Cm = ^ f ^l^™] ^ 1 and hence to 

I |2. _ ( SL[Rm] y^^ 2.^ 

The Spectral Renormalization method is therefore given by the iterations 

In this work, we use the SRM to solve l\17\i for the cases d = 1,2,3 without imposing radial 
symmetry. Alternatively, one might have solved the radial equation (120p using a modified Hankel- 
like transform instead of the Fourier Transform. Our main reason for not doing so is the convenience 
and cost-effectiveness of using the Fast Fourier Transform. We also note that our non-radially- 
symmetric method produced a radially-symmetric solution, which suggests that the ground state 
is radially symmetric. 
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A Compactness Lemma 

Here we provide an extension of the Compactness Lemma for -f^^adiai functions [32j, to the case 

Lemma 13. (Compactness Lemma) Let d> 2 and let a > be in the H'^-suhcritical regime ([6]). 
Then, the embedding Hl^^^^^{W^) ^ L'^^''+'^\W^) is compact, i.e., every bounded sequence Un' £ 
-^radial (^'^) ^ subscquencc Un which converges strongly in L'^^'^^^\W^) . 

Proof. If ||iin'|l_ff2 < M, then the sequence Un' has a subsequence n„ which converges weakly to u 
in H^. Since the limit of radial functions is a radial function, u G ^^radiai- addition, since for 
any bounded domain Q, the embedding H'^(Q) L'^(Q) is compact, there is a subsequence which 
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converges strongly to u in L'^(Q), i.e., lim„_>oo /q \un — n|^(ix = 0. From the Gagliardo-Nirenberg 
inequality on the bounded domain Q, see fl5 \ 116 1 [26 ] . 

and since 1 > cr^^ in the -ff^-sub critical case, it follows that Un ^ u strongly in L^*^'^^^^ ($7), so 
that 

lim / \un — n|^'-'^^^^ dx = 0. 



n— >oo 



n 



Next, Strauss radial Lemma |32j for functions gives that Vpe > 1 and n, 



, , {d-l)(T' 



SO that Ve3pe s.t. Vn 

/ |n„p('^+i)dx<e. 

.7 \x\>pe 

Finally, since 



□ 



\\Un - n||^2{<T+l)(]Rd) < \\Un - ■w|Il2{<t+i)(|x|<p^) + ll^n || l2(<t+i) (|x|>p^) + II'"IIl2{-+1){|x|>p,) 

the convergence in R'^ is obtained. 

B Proof of Corollary [6] 

Proof. From the proof of Theorem U 

and since limfc^oo o,{tk) / K'^k) = c«, we have that VAf > 

liminf ll0fc(?')lli2(r<Af) < 1™™^ ll'^fcWllL2{r<a(t,)/i{tfc)) ' 
/c — ^ oo ^ oo 

Since (ii.fr)^^', it follows that (bhir) ^, and so 

ll^llL2(r<Af) < l™inf ll'/'fcW||i2(^<M) 

This is true VM, and so 

Per < < l^i^f IIV'(ifc,0lli2(.<a(t,)) • 

For the second result, since ||^'||/^2 > -Per; it follows that for all e > there exist K > such 
that 

ll^llL2(.>i^) > (l-e)^'cr. 
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Therefore, since 

Il^(ifc,0lli2(r<7^./(t,)) = \\Mr)\\l2(^r<K) ' 

a similar argument as in the previous section gives 

(l-e)^'cr < \\'^\\h(r<K) < l™inf ||V'(ife,0llL2(r</^.Z(t,)) • 

□ 



C Gagliardo Nirenberg inequality for functions 



In L^-critical case ad = 4, which is always in the if^-subcritical regime 
Gagliardo-Nirenberg inequality in is [XB } [TUl 

Lemma 14. (Gagliardo-Nirenberg inequality) Let ad = 4, and let f G H'^(R'^ 

|^[^+;;<i?.,,llA/||^||/||2^ 



the appropriate 



then 



(49) 



We note that the ground-state R of equation ^ is the minimizer of the Gagliardo-Nirenberg 
inequality p!2], and that its norm, the critical power, satisfies 



\R\\ 



a + l 

BcrA 



1/a 



Hence, the Gagliardo-Nirenberg inequality implies the following Corollary: 
Corollary 15. Let f G and ad = 4, then 



H[f]> 



1 



SO that 



< P 

2 — cr 



l|A/||^ 
H[f] > 0. 



(50) 



D WKB analysis of eq. ( Md) . 



As p — > oo, the nonlinear term in (I44cp becomes negligible, and (I44cp reduces to 



B{p) - ApB + ih-^ -B + pB^ 



a 



0, 



(51) 



where 



(d-l)(d-3) (d-l)((i-3) 2 , 2(d-l) 3 4 

3 <^p i- 2 '^p "T ^ P P 
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In order to apply the WKB method, we substitute B{p) = exp{w{p)), and expand 

w{p) ~ woip) +'wi{p) + ... . 

Substituting wq{p) = apP and balancing terms shows that p = 4/3, and that the equation for the 
leading-order, the O {p'^^^) terms, is 

3 



KV'=(f) P = ^»', 



Therefore, 



3, i2L+i^2fc 3 j . VS + i -Vs + i 



The equation for the next order, the 0(1) terms, is 

1 + 2dib^ = zftS + 3pw[^ , 

implying that 

w, = l(^^{l-ad)-^yogp. 

The next-order terms are can be neglected. We therefore obtain the three 

solutions 

B2{p) ~ -^^^exp(^-i^6//3-i^log(p)^, (52a) 

Bsip) - ;^(,,„,) ^ exp(+.-5//3-.-log(p)j, (52b) 

exp (- 3^6p4/3) / 3 , 1 \ 
Bdp) ^ exp (+^-6//3 _ j . (52e) 

Since (j5ip is a fourth order ODE, another solution is required. To obtain the fourth solution, we 
substitute wq ~ /31og(p) in ([?T|) and obtain that the equation for the leading-order, the 0(1) terms, 
is 

+ (- + /?) =0, 



and that the next order is o(l) and can be neglected. The fourth solution is therefore 

Blip) ^ p'-^'l^. (52d) 
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E Application of Noether Theorem for the BNLS 

The Lagrangian density of the BNLS is 

Z i + <T 

We cite here Noether's Theorem, as given in |33j : 



2{a+l) 



Theorem 16 (Noether Theorem). If the action integral JJ Cdxdt is invariant under the infinites- 
imal transformation 

t h-f t = t + 5t(x, t, -i/;), 
X ^ X = X + Jx(x, i, -0), 
ijj I— > ^p = ^p + 6ip{x,t,'ip), 



then 



dC dC 

— iiJtSt + • 5x - #) + TT-r 0^*1 St + V^* -dx-sr)- C5t 



dx 



(53) 



is a conserved quantity. 



For example, the BNLS action integral is invariant under the phase-multiplication ■0(i, x) 
e*^^/'(t,x). In this case, 5t = 0, 5x = 0, 8ij) = iip, and so Theorem [16] implies that the integral 



dC dC 

— iiJtSt + Vi,.Sx-Si;) + — {rt5t + Vr ■ <5x - <5V*) - 
diJt oipl 

^V* (V'i • + • - iV') - ^t/j {ipt ■ + VV'* • + iV*) - C-Q 

2 

2 ' 



dx 
dx 



i.e., the power, is a conserved quantity. Other conservation laws can be found in a similar manner. 
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